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ABSTRACT 

Several  recent  studies  of  internal  gravity  waves  have  been  expressed  in  a  Lagrangian  reference  frame,  motivated 
by  the  observation  that  in  this  frame  the  dispersion  relation  then  excludes  the  Eulerian  Doppler  shifting  term 
due  to  a  background  flow.  Here  the  dispersion  relation  in  a  Lagrangian  reference  frame  is  explicitly  derived  for 
a  background  flow  and  background  density  that  are  slowly  varying  with  respect  to  the  waves,  but  are  otherwise 
arbitrary  functions  of  space  and  time.  Two  derivations  are  given,  both  yielding  the  same  result.  The  first  derivation 
involves  a  transformation  of  the  dispersion  relation  from  Eulerian  to  Lagrangian  coordinates,  while  the  second 
derivation  involves  a  wave-packet  analysis  of  the  equations  of  motion  directly  in  Lagrangian  coordinates.  The 
authors  show  that,  although  the  Eulerian  Doppler  shifting  term  is  removed  from  the  dispersion  relation  by  the 
transformation  of  the  frequency  when  passing  from  an  Eulerian  to  a  Lagrangian  reference  frame,  a  dependence 
on  the  background  shear  is  then  introduced  by  the  transformation  of  the  wavenumber.  This  dependence  on  the 
background  shear  is  the  term  that  accounts  for  wave  refraction  in  the  Lagrangian  frame,  and  its  role  has  apparently 
not  been  fully  appreciated  in  the  aforementioned  previous  studies. 


1.  Introduction 

There  has  been  a  recent  revival  of  interest  in  for¬ 
mulating  a  theory  for  internal  gravity  waves  in  a  La¬ 
grangian  frame  of  reference,  rather  than  the  usual  Eu¬ 
lerian  frame  of  reference.  Eollowing  the  work  of  Allen 
and  Joseph  (1989),  Hines  (2001),  and  Chunchuzov 
(2002)  have  argued  that  the  nonlinear  terms  in  the  La¬ 
grangian  form  of  the  equations  are,  in  some  cases,  of 
less  importance  than  the  corresponding  nonlinear  terms 
in  the  Eulerian  form  of  the  equations.  In  effect,  the 
implication  is  that  simply  using  a  coordinate  transfor¬ 
mation  from  the  Eulerian  frame  to  the  Lagrangian  frame 
can  remove  some  of  the  effects  of  nonlinearity.  This  is 
an  interesting,  albeit  provocative,  suggestion.  It  is  of 
course  true  that,  in  some  special  cases,  a  system  of 
partial  differential  equations  can  be  totally  or  partially 
explicitly  linearized  by  a  coordinate  transformation,  but 
this  is  not  so  in  general  for  the  transformation  from 
Eulerian  to  Lagrangian  coordinates. 
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The  above  studies  are  concerned  with  an  internal 
wave  field  in  a  background  composed  of  other  internal 
waves.  In  each  of  these  studies,  the  Lagrangian  internal 
wave  dispersion  relation  is  taken  to  have  the  same  func¬ 
tional  form  as  the  well-known  Eulerian  internal  wave 
dispersion  relation,  but  without  the  Doppler  shifting 
term  [see,  for  instance,  (2.3)  of  Allen  and  Joseph  (1989) 
or  (59)  below].  Hence,  such  a  Lagrangian  dispersion 
relation  loses  a  crucial  dependence  on  the  background 
flow  and,  in  particular,  cannot  account  for  refraction  by 
the  background  flow.  Thus,  in  particular,  Hines  (2002) 
has  suggested  that  under  appropriate  conditions  (see  sec¬ 
tion  5),  there  are  consequent  simplifications  for  the  ray 
tracing  of  internal  waves  if  the  rays  are  expressed  in 
Lagrangian  coordinates. 

These  issues  have  motivated  us  to  reexamine  the 
structure  of  the  internal  wave  dispersion  relation  in  a 
Lagrangian  reference  frame  vis-a-vis  the  corresponding 
internal  wave  dispersion  relation  in  an  Eulerian  refer¬ 
ence  frame.  Hence,  in  this  paper,  we  explicitly  and  sys¬ 
tematically  derive  the  Lagrangian  dispersion  relation  for 
small-amplitude  internal  waves  in  the  presence  of  a  gen¬ 
eral  background  flow  of  finite  amplitude.  Lor  simplicity, 
we  consider  two-dimensional  flow  and  ignore  the  effects 
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of  the  earth’s  rotation.  Two  derivations  are  given,  which 
lead  to  the  same  result.  The  first  involves  a  transfor¬ 
mation  of  the  dispersion  relation  from  Eulerian  to  La- 
grangian  coordinates  (section  3a),  and  the  second  in¬ 
volves  a  wave-packet  asymptotic  analysis  of  the  La- 
grangian  equations  of  motion  (section  3b).  A  general, 
slowly  varying,  space-  and  time-dependent  background 
is  treated  hrst.  Then,  for  illustration,  we  specialize  to  a 
horizontal  background  flow  that  varies  only  in  the  ver¬ 
tical  direction  (section  4).  The  main  conclusion  from 
our  analysis  is  that  there  is  a  complete  equivalence  be¬ 
tween  the  Eulerian  and  Lagrangian  wave  dispersion  re¬ 
lations  and  that  internal  wave  rays  that  refract  in  an 
Eulerian  frame  will  refract  in  a  Lagrangian  frame,  and 
vice  versa. 

2.  The  Eulerian  dispersion  relation 

The  Eulerian  coordinates  are  x  =  (x,  z)  with  z  vertical. 
We  use  an  overbar  to  denote  quantities  associated  with 
the  background.  Eor  example,  the  background  velocity 
u(x,  t)  =  (u,  W)  and  the  background  density  p(x,  f). 
The  background  flow  is  incompressible  and  satishes  the 
equations  for  conservation  of  density,  but  we  allow  for 
the  possible  presence  of  arbitrary  body  forces  in  the 
momentum  equations  to  maintain  the  background  flow. 

Linearizing  about  this  background,  the  equations  for  the 
associated  perturbations  u,  w,  p,  and  pressure  p  are 


p{u,  +  uu^  -f  wuj  +  +  ■  ■  ■  =  0,  (1) 

p(w,  +  uw^  +  wwj  +  p.  +  gp  +  ■  ■  ■  =  0,  (2) 

p,  +  Mp,  +  vvpj  +  up^  -I-  wp;  =  0,  and  (3) 

+  Wj  =  0.  (4) 

Here,  subscripts  denote  partial  derivatives,  and  g  is  the 
acceleration  due  to  gravity.  In  the  momentum  equations 
{ •  •  • }  denote  terms  that  can  be  neglected  at  leading  order 
since  they  involve  derivatives  of  u,  w.  More  formally, 
we  could  introduce  the  slow  variables  X  =  ex,  Z  =  ez, 
T  =  et,  where  e  is  a  small  parameter  characterizing  the 
slow  variation  of  the  background  with  respect  to  the 
wave  held,  so  that  u  =  u{X,  Z,  T),  etc.  The  derivatives 
are  then  0(e)  with  respect  to  the  wave  frequency.  But 
it  is  important  here  to  note  that  g'pj'p  and  gpjp  are 
(9(1)  quantities  with  respect  to  the  square  of  the  wave 
frequency.  The  reason  for  this  is  that  for  these  internal 
waves  the  wave  frequency  scales  with  (g/H)"^,  where 
//  is  a  scale  height  for  the  density  stratihcation.  To  avoid 
excessive  notation,  we  will  not  formally  introduce  these 
slow  variables  into  the  derivation,  although  their  pres¬ 
ence  is  always  understood. 

Introducing  the  perturbation  streamfunction 

M  =  -i//,,  W  =  lA,,  (5) 

and  eliminating  the  pressure  from  the  momentum  equa¬ 
tions  gives,  at  leading  order. 


p[(V^<A),  +  +  vv(V^^)J  +  gp^+  ...  =0, 

(6) 

p,  +  +  wPz  -  'AjA,  +  4>.Pz  =  0. 

(7) 


We  let 


ip,  p  a(x,  -I-  C.C.,  (8) 

where  a  is  the  amplitude,  6  is  the  phase  for  a  wave 
packet,  and  c.c.  denotes  complex  conjugate.  More  for¬ 
mally,  a  =  a  (X,  Z,  T)  is  slowly  varying  on  the  same 
scale  as  the  background,  whereas  the  phase  6  =  8(X, 
Z,  T)/e  is  rapidly  varying.  The  Eulerian  wavenumber 
and  frequency  are  then  dehned  by 

k  =  (k,  m)  =  V0,  (9) 

w  =  -e,.  (10) 

Because  the  phase  is  rapidly  varying  with  respect  to  the 
background,  at  leading  order  the  derivatives  of  the  held 
variables  are  taken  with  respect  to  the  phase  only.  Sub¬ 
stituting  (8)  into  (6)-(7)  then  leads  to  the  Eulerian  dis¬ 
persion  relation 

^2  =  (k^N^  -  kmM^)l(k^  -t-  m"),  (11) 

with 


w  =  w  —  ku  —  mw.  (12) 

Here  =  —gji-jj}  and  =  —g'pj'p.  Often  the  hor¬ 
izontal  density  gradient  is  ignored;  that  is,  M  =  0. 

An  equivalent  representation  for  the  Eulerian  disper¬ 
sion  relation,  and  one  that  is  more  convenient  for  our 
derivation  of  the  Lagrangian  dispersion  relation  in  sec¬ 
tion  3a,  is  obtained  by  using  the  fact  that  the  background 
density  p  satishes  the  equation  for  the  conservation  of 
density.  This  can  be  solved  to  give  p  =  pjz  —  ^).  Here 
Po(z)  is  the  undisturbed  density  held  (in  the  absence  of 
any  how)  and  ^  is  the  vertical  displacement  induced  by 
the  background  how,  which  satishes  the  equation: 

(13) 

We  then  use  the  relations 

Pz  =  Pod  -  P.  =  Po(-?d,  (14) 

where  the  prime  indicate^differentiation  of  p„  with  re¬ 
spect  to  its  argument  (z  —  ^).  Consequently,  the  Eulerian 
dispersion  relation  can  be  rewritten  as 


N^kk/(k^  +  mj. 

(15) 

k  +  mj^  -  kj^. 

(16) 

-gp'JPo. 

(17) 

Note  that  N  =  N(z  —  O  so  varies  on  the  same 
temporal  and  spatial  scales  as  the  background  how. 
The  Eulerian  ray  equations  (Lighthill  1978)  are  read- 
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ily  obtained  from  the  Eulerian  dispersion  relation  w  = 
(u(k,  X,  f)  in  (15)  and  are  given  by 


dx 

dt 


=  Vjw, 


(18) 


Here  didt  is  the  rate  of  change  following  the  ray  at  the 
local  group  velocity  =  (wj,,  w„,).  Along  an  Eulerian 
ray,  we  also  have  dwidt  =  u>,. 


3.  The  Lagrangian  dispersion  relation 


In  Lagrangian  coordinates  x'  =  (x’,  z'),  the  equations 
of  motion  are  (e.g.,  Lamb  1932) 


Po(v„x,.  +  +  gz,.)  +  P,.  =  0,  (19) 

Po(v„x-,  +  z„z,.  +  gz,.)  +  p,.  =  0,  (20) 


P  =  Poiz') 


d(x,  z) 
d{x',  z’) 


=  1. 


(21) 


To  derive  the  Lagrangian  dispersion  relation,  we  can 
apply  an  analogous  wave-packet  analysis  to  these  La¬ 
grangian  equations,  or  we  can  transform  the  Eulerian 
dispersion  relation  to  the  Lagrangian  frame  of  reference. 
We  give  both  derivations,  starting  with  the  latter. 


a.  Derivation  by  transformation  of  the  Eulerian 
dispersion  relation 

This  derivation  is  based  on  the  mapping  between  Eu¬ 
lerian  and  Lagrangian  coordinates,  whose  Jacobian  is 
given  in  (21).  The  Lagrangian  momentum  equations 
(19)-(20)  are  not  explicitly  used  in  this  derivatiom  _ 
We  introduce  the  particle  displacements 
such  that 

X  =  x'  +  ^(x’,  z',  t)  -t-  ^(x',  z' ,  t),  and  (22) 

z  =  z'  +  Cix',  z',  f)  +  I{x’,  z',  t).  (23) 

Ir  i^  sufficient  here  to  consider  only  the  displacements 
f  due  to  the  background  flow.  To  retain  the  pertur¬ 
bation  displacements  ^  would  invalidate  the  lineari¬ 
zation  procedure  as  applied  in  both  the  Eulerian  and 
Lagrangian  frames.  In  particular,  their  retention  would 
lead  to  their  occurrence  in  the  Lagrangian  dispersion 
relation,  clearly  inconsistent  with  a  linear  theory.  Lor 
the  rest  of  this  subsection,  we  therefore  neglect  f  L  in 
(22)-(23). 

Lor  the  background  flow,  the  particle  displacement 
helds  are  related  to  the  Eulerian  background  velocity 
held  through  the  equations 

1,  =  u{x’  -f  z'  +  I,  i),  and  (24) 

=  vv(x'  +  I,  z'  +  ?,  0,  (25) 

with  ^  =  0,  ^  =  0  at  f  =  0.  Lormally,  the  background 
displacements  depend  on  the  slow  variables  X'  =  ex', 
Z'  =  ez’  ,T  =  et,  which  are  the  Lagrangian  counterparts 
of  the  Eulerian  slow  variables  X  =  ex,  Z  =  ez,  T  =  et 


dehned  in  section  2.  Note  that  the  displacements  are 
formally  of  0(1/ e),  although  of  course  their  hrst  deriv¬ 
atives  are  0(1).  As  in  the  Eulerian  case,  we  will  not 
explicitly  use  these  slow  variables  here,  but  their  pres¬ 
ence  is  always  understood. 

It  will  be  useful  for  the  present  derivation  to  construct 
the  matrix  B,  dehned  by 


B  = 


1  + 


(26) 


Incompressibility  of  the  background  how  implies  that 
detB  =  1. 

To  transform  the  Eulerian  dispersion  relation  into  La¬ 
grangian  coordinates,  we  must  transform  not  only  k,  m. 
Cl)  but  also  the  Eulerian  derivatives  that  occur  in 

the  dehnition  of  k  in  (16).  Starting  with  this  latter  trans¬ 
formation  we  hrst  note  that 


so  that 

x;  =  Z,.,  z'  =  x,,.,  x'  =  -X,.,  z^  =  -z,..  (29) 

These  four  equations,  on  using  (22)-(23)  for  the  back¬ 
ground  how  only,  imply  that 

-I,  =  J,.,  -J,  =  =  I,., 

L  =  L-  (30) 

Hence  k  in  (16)  becomes 

k  =  k  +  m|,.,  -f  kf  ^..  (31) 

Next  we  transform  the  Eulerian  wavenumber  and  fre¬ 
quency  k.  Cl).  The  Lagrangian  wavenumber  and  fre¬ 
quency  will  be  denoted  by  k',  co'.  We  start  with  the 
phase  9  in  (8),  which  transforms  from  Eulerian  to  La¬ 
grangian  coordinates  according  to 

9(x,  t)  =  0[x(x',  t),  t]  =  0'(x',  t).  (32) 

Dehning  V'  =  (5,.,,  we  have 

k'  =  V'0'  Cl)'  =  -9',.  (33) 

Using  =  0,.x,.  -t-  fl.Zt-  and  a  similar  expression  for 
9'.,,  we  arrive  at  the  wavenumber  transformation 

k'  =  kx„.  -t-  mZ:,',  (34) 

m'  =  kx^.  +  mz^’.  (35) 

Equivalently, 

k'  =  BTk,  (36) 

where  superscript  T  indicates  the  transpose.  Using  the 
second  expression  for  B  in  (26),  we  have 
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k'  =  k  +  -f 

(37) 

m'  =  m  +  kf^.  m^^,. 

(38) 

Lrom  (31)  and  (37)  we  now  see  that  k'  =  k;  that 
is  in  fact  the  Lagrangian  horizontal  wavenumber. 
The  inverse  relations  are 

is,  k 

k  =  k'  +  k'i^,  - 

(39) 

m  =  m’  —  -I- 

(40) 

Lor  the  frequency. 

co'  =  -e,  -  e,x,  -  0^z„ 

(41) 

=  0)  —  ku  —  mw. 

(42) 

=  d). 

(43) 

Thus  the  Eulerian  dispersion  relation  (15)  can  be 
written  as 


0)2  =  mkk’lik^-  +  m2).  (44) 

Substituting  for  w,  k,  m  from  (39)-(43)  gives  the  La- 
grangian  dispersion  relation 

^  _  mk'(k^  +  k~C,  -  m'J^) 

(k'  +  k'C^.  -  +  (m’  -  k'^,  + 

(45) 

The  derivatives  of  ^  appearing  here  are,  by  (24)-(25), 
functions  of  the  background  shear. 

The  Lagrangian  ray  equations  are  the  counterpart  of 
the  Eulerian  ray  equations  (18).  That  is,  given  the  La¬ 
grangian  dispersion  relation  o)'  =  o)'(k',  x',  t)  in  (45), 
the  Lagrangian  ray  equations  are 


Po[e,(i  +  L')  +  LL'  +  gL']  +  P.'  =  0,  (48) 

+  Ui  +  ?.■)  +  gU  +  P,  =  0,  (49) 

^,,(1  +  ?,,)  +  ^,,(1  + 1,0  -  -  CA,  =  0.  (50) 

Eliminating  the  pressure  from  (48)-(49)  gives,  to  lead¬ 
ing  order, 

iL'  -  Un  +  LA.'  +  LA.’ 

-  LAL'  -  LA.'  +  =  0.  (51) 

We  seek  solutions  of  the  form 

^  oc  a'{x',  Oe'®'**'’')  -t-  C.C.,  (52) 

with  k'  =  V'P'  and  w'  =  —6',  as  in  (33).  Substitution 
into  (50)-(51)  gives 

w'2(m^  -  kO  +  N^-k'C  =  0,  (53) 

mC  +  k^=0,  (54) 

where  we  have  used  the  wavenumber  relation  (36)  to 
write  the  result  in  a  compact  form.  Eliminating  ^  or  ^ 
from  these  equations,  we  obtain 

w'2  =  N^kk'Hk^  +  m2),  (55) 

which  agrees  with  (43)  and  (44). 

4.  A  steady  horizontal  background  flow 

We  now  apply  the  general  Lagrangian  dispersion  re¬ 
lation  just  obtained  to  the  special  case  in  which  the  back¬ 
ground  flow  M  is  horizontal  and  varies  with  depth  only. 
The  background  displacements  are  then  given  by 

I  =  u{z')t,  1  =  0.  (56) 


Along  a  Lagrangian  ray,  du>' Idt  =  u>’,.  Using  the  trans¬ 
formations  (22)-(23),  for  the  background  flow  alone, 
along  with  (36)  and  (42)-(43),  it  can  now  be  shown  that 
there  is  a  complete  equivalence  between  these  Lagrang¬ 
ian  ray  equations  (46)  and  the  Eulerian  ray  equations 
(18).  The  demonstration  is  straightforward,  but  note  the 
relation 


B,  = 


(47) 


obtained  by  differentiation  of  (24)-(25)  with  respect  to 


b.  Derivation  by  an  asymptotic  analysis  of  the 
Lagrangian  equations 

This  derivation  first  requires  a  linearization  of  the 
Lagrangian  equations  (19)-(21)  about  the  background 
flow,  using  the  background  and  perturbation  displace¬ 
ments  defined  in  (22)-(23).  The  result  is 


Lrom  (34)-(35),  the  relation  between  the  Eulerian  and 
Lagrangian  wavenumbers  is  then 

k  =  k',  m  =  m'  —  k'u^.t.  (57) 

The  Lagrangian  dispersion  relation  (45)  reduces  to 

w'2  =  k’^N^lk'^  +  (m'  -  k'u^.tf].  (58) 

The  background  shear  is  small  compared  with 
oi'since,  when  expressed  in  terms  of  the  slow  variables, 
it  is  ed^..  This  does  not  imply,  however,  that  we  can 
neglect  k'u^.t  in  comparison  with  m'  in  (58).  Indeed,  in 
terms  of  the  slow  variables  this  term  is  k'u^.T  and  is 
0(1)  with  respect  to  the  small  parameter  e.  Alternatively 
we  could  estimate  k'u^.t  as  k'ulc^.  Here  is  the  vertical 
component  of  group  velocity  (which  for  this  problem 
has  the  same  value  in  both  the  Eulerian  and  Lagrangian 
frames).  We  have  used  the  estimates  that  along  a  ray 
path  u,,z'  ~  u  and  z' It  ~  c^.  Thus,  k’u^.t  is  an  0(1) 
quantity  and  of  the  same  order  as  m'. 

It  is  clear  from  the  study  of  Hartman  (1975)  that  the 
term  k’u^t  in  (58)  can  be  important  for  the  description 
of  rays  in  the  Lagrangian  frame.  He  considers  the  par¬ 
ticular  case  of  constant  N  and  constant  background  shear 
Then  m',  as  well  as  k',  are  constant  along  a  La- 
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grangian  ray,  but  w'  is  not  constant.  This  follows  from 
the  Lagrangian  dispersion  relation  (58),  which  depends 
explicitly  on  t  but  not  on  z'  when  m..  is  constant.  For 
this  particular  case,  the  Lagrangian  rays  are  straight  lines 
in  (x',  z'),  but  the  group  velocity  changes  along  the  ray 
due  to  the  change  in  w'.  The  Lagrangian  ray  reflects 
from  the  same  turning  point  height  (where  co'  N) 
and  asymptotes  toward  the  same  critical  layer  height 
(where  w'  0)  as  for  the  corresponding  Eulerian  ray. 

Note  that  the  coordinates  used  by  Hartman  (1975)  are 
Lagrangian  in  the  sense  that  we  have  defined  here.  Hart¬ 
man’s  coordinates  are  x'  =  x  —  u(z)t  and  z'  =  z,  that 
is,  the  position  that  moves  with  the  velocity  of  the  back¬ 
ground  flow.  As  the  notation  indicates,  these  are  the 
same  as_our  Lagrangian  coordinates  x'  =  x  —  f  and  z' 
=  z  —  ^  when  our  coordinates  are  specialized  to  this 
case  using  (56). 

5.  Discussion 

The  main  result  of  this  paper  is  the  Lagrangian  dis¬ 
persion  relation  (45)  for  internal  gravity  waves.  This 
differs  from  the  form  used,  for  instance,  by  Allen  and 
Joseph  (1989),  which  in  present  notation  is  (apart  from 
the  Coriolis  term  and  ignoring  the  third  dimension) 

w'^  =  N^k'y(k'^  +  m'^).  (59) 

Unlike  this  form,  which  is  appropriate  for  a  uniform 
background  without  any  shear,  our  Lagrangian  disper¬ 
sion  relation  (45)  depends  on  the  background  straining 
held.  That  is,  it  contains  terms  involving  the  derivatives 
of  the  particle  displacements  ^  and  ^  associated  with 
the  background  flow.  These  terms  are  functions  of  the 
background  shear  through  (24)-(25),  and  their  occur¬ 
rence  in  the  Lagrangian  dispersion  relation  then  ac¬ 
counts  for  wave  refraction  in  the  Lagrangian  frame. 

Allen  and  Joseph  (1989)  used  Lagrangian  coordinates 
to  describe  a  spectrum  of  internal  waves,  with  each  wave 
component  of  the  Lagrangian  spectrum  satisfying  (59). 
When  this  Lagrangian  spectrum  is  transformed  to  the 
Eulerian  frame,  the  resulting  high-wavenumber  com¬ 
ponents  of  the  Eulerian  spectrum  were  found  to  be  non¬ 
wavelike.  That  is,  these  components  did  not  satisfy  the 
Eulerian  dispersion  relation  (15).  Our  work  suggests  a 
possible  factor  contributing  to  this  result,  in  that  Allen 
and  Joseph’s  Lagrangian  dispersion  relation  (59)  and 
the  Eulerian  dispersion  relation  (15)  are  not  equivalent 
to  each  other  under  the  Lagrangian-to-Eulerian  coor¬ 
dinate  mapping  except  in  the  limit  of  no  background 
shear.  Eor  these  high  wavenumbers,  this  limit  is  not 
appropriate  because  the  high  wavenumbers  are  sub¬ 
jected  to  strong  background  shears  resulting  from  the 
longer  waves.  A  comprehensive  discussion  of  this  issue 
of  the  high-wavenumber  Eulerian  spectrum  can  be 
found  in  the  recent  review  by  Eritts  and  Alexander 
(2003). 

Hines  (2002)  examines  another  issue.  He  considers  a 
“test  wave”  in  a  background  of  waves  and  asks  the 


question:  when  can  the  interaction  of  the  test  wave  with 
the  background  be  ignored?  Hines  then  finds  that  the 
criteria  for  neglecting  the  interaction  with  the  back¬ 
ground  depends  on  the  test-wave  wavenumber  in  dif¬ 
ferent  ways  according  to  whether  one  uses  an  Eulerian 
or  a  Lagrangian  coordinate  system  (see  Eig.  1  of  Hines 
2002). 

In  the  Lagrangian  frame,  the  Hines  criteria  for  neg¬ 
ligible  interaction  with  the  background  involve  his 
quantities  S^  i,  533^ which  correspond  respec¬ 
tively  to  our  quantities  The  conditions 

on^^i  1,  in  Hints’  (3.17),  (3.22),  correspond  here  to 
I  I  ^1  and  I  I  <C  I  m' Ik’  I .  The  conditions  on 
53 1,  533  in  Hines’ J^3. 21),  (3.18),  correspond  here  to  1^^.  I 
<C  \k’/m'  I  and  \^^,\  <C  1.  Under  these  conditions,  the 
background  terms  in  our  Lagrangian  dispersion  relation 
(45)  become  small,  so  in  this  respect  our  work  is  con¬ 
sistent  with  that  of  Hines  (2002). 

In  an  Eulerian  frame  of  reference,  the  corresponding 
conditions  under  which  the  background  flow  can  be  ig¬ 
nored  are  readily  obtained  from  (15).  Thus,  we  now 
require  the  conditions  under  which  w  ~  w,  or  equiva¬ 
lently  I  k  •  u  I  <C  w,  and  under  which  k  ~  k.  On  using 
the  definitions  (16)  and  (15),  these  can  be  expressed^in 
the  form  lul  <CNIkl/lkU,  1^,1  <C  \klm\  and  I  ^ J 
<C  1 .  Again,  these  conditions  are  consistent  with  those 
obtained  by  Hines  (2002)  for  the  neglect  of  the  effects 
of  the  background  in  an  Eulerian  reference  frame.  How¬ 
ever,  we  point  out  here  that  the  criterion  I  k  •  u  I  <C  <it> 
essentially  estimates  advection  by  the  background  flow, 
rather  than  refraction,  per  se.  To  take  account  of  re¬ 
fraction,  one  needs  to  estimate  the  variation  of  k  ■  u 
along  a  ray.  When  this  is  done,  the  criterion  is  replaced 
by  I  5u  I  ^  M I  k  I ,  where  Su  is  a  measure  of  the  spatial 
variability  in  u.  Eurther,  we  also  need  to  ensure  that  A(z 
—  C)  does  not  vary  significantly  _with  [.  This  can  be 
achieved  by  requiring  either  that  ^1  and  I  ^  ^ 

1 ,  or  that  N  itself  is  approximately  constant. 

Hines  (2002)  argues  that  there  is  an  advantage  in 
using  a  Lagrangian  frame  when,  as  judged  by  the  above 
criteria,  the  interaction  between  the  Lagrangian  ray  and 
the  background  is  weak  (see  appendix  C  of  Hines  2002). 
But  note  that  the  above  criteria  are  local  and  wavenum¬ 
ber  dependent.  They  can  be  used  to  describe  the  con¬ 
ditions  under  which  the  refraction  is  locally  minimal, 
either  in  an  Eulerian  frame  on  the  one  hand  or  in  a 
Lagrangian  frame  on  the  other  hand.  However,  the  pres¬ 
ence  of  variations  in  the  background,  no  matter  how 
weak,  will  cause  some  refraction  and,  since  the  criteria 
for  minimal  refraction  are  wavenumber  dependent,  it  is 
not  clear  whether  the  refraction  can  necessarily  remain 
minimal  along  a  ray  over  the  time  scales  of  interest. 
Indeed  Hines  (2002,  p.  26)  notes  that  such  conditions 
for  minimal  refraction  may  be  breached  locally  and  need 
to  be  continually  confirmed.  On  the  other  hand,  our 
present  formulation  based  on  the  full  Lagrangian  dis¬ 
persion  relation  (45)  is  not  subject  to  any  such  restric¬ 
tion.  Numerical  ray  tracing  based  on  the  full  Lagrangian 
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dispersion  relation  (45)  might  help  to  resolve  this  ques¬ 
tion,  and  to  test  the  criteria  derived  by  Hines  (2002)  and 
described  above. 
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